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Abstract 

We present detailed results for the model-independent low momentum nucleon- 
nucleon interaction V\ ow k- By introducing a cutoff in momentum space, we separate 
the Hilbert space into a low momentum and a high momentum part. The renormal- 
ization group is used to construct the effective interaction v\ ow k in the low momen- 
tum space, starting from various high precision potential models commonly used in 
nuclear many-body calculations. With a cutoff in the range of A ~ 2.1 fm _1 , the new 
potential V\ ow k is independent of the input model, and reproduces the experimental 
phase shift data for corresponding laboratory energies below E^ ~ 350 MeV, as 
well as the deuteron binding energy with similar accuracy as the realistic input po- 
tentials. The model independence of V\ ow k demonstrates that the physics of nucleons 
interacting at low momenta does not depend on details of the high momentum dy- 
namics assumed in conventional potential models. V\ ow k does not have momentum 
components larger than the cutoff, and as a consequence is considerably softer than 
the high precision potentials. Therefore, when V\ ow k is used as microscopic input in 
the many-body problem, the high momentum effects in the particle-particle chan- 
nel do not have to be addressed by performing a Brueckner ladder resummation or 
short-range correlation methods. By varying the cutoff, we study how the model 
independence of Vj ow k is reached in different partial waves. This provides numerical 
evidence for the separation of scales in the nuclear problem, and physical insight 
into the nature of the low momentum interaction. 
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1 Introduction 



A major challenge of nuclear physics lies in the description of finite nuclei 
and nuclear matter on the basis of a microscopic theory. By microscopic, it 
is understood that the input nuclear force is, as best as possible, based on 
the available free-space data, e.g., the elastic nucleon-nucleon and nucleon- 
deuteron scattering phase shifts. The theoretical predictions are then obtained 
using a systematic many-body approach, which is reliable for strongly inter- 
acting systems. 

At low energies and densities, one typically starts by treating protons and 
neutrons as nonrelativistic point-like fermions. However, as opposed to the 
Coulomb interaction in electronic systems, the nucleon-nucleon interaction re- 
mains unknown from the underlying theory of the strong interactions, Quan- 
tum Chromodynamics (QCD). Thus, phenomenological meson-exchange mod- 
els of the two-nucleon force are commonly used as input in many-body calcu- 
lations. Moreover, three- (and possibly four-) body forces are needed to obtain 
accurate results for few-body systems. Recent Green's function Monte Carlo 
calculations using a realistic two-body interaction, the Argonne fig potential, 
have shown that three-body forces typically contribute a net attraction of 
~ 20% to the low-lying states of light nuclei, A ^ 10 [1,2]. 

However, already at the lowest level of the two-body interaction, there are sev- 
eral quite different potential models commonly used. Therefore, it is a natural 
first step to ask whether it is possible to remove or to minimize this model 
dependence. In this work, we expand on previous results [3,4] and demonstrate 
that this is indeed possible, as long as one restricts the interaction to momenta 
which are constrained by the experimental scattering data, and explicitly in- 
cludes the long-range part of the interaction. 

There are a number of high precision models of the nucleon-nucleon force Vnn, 
which we will refer to as the "bare" interactions. All realistic potentials are 
given by the one-pion exchange (OPE) interaction at long distances, but vary 
substantially in their treatment of the intermediate-range attraction and the 
short-range repulsion. We briefly summarize the different frameworks used to 
generate the intermediate and short-range parts of the realistic potential mod- 
els [5,6]. The Paris potential [7] calculates the intermediate-range contributions 
to the nucleon-nucleon scattering amplitude from two-pion exchange using dis- 
persion theory, and then assumes a local potential representation consisting of 
several static Yukawa functions. The cu-meson exchange at short distances is 
included as part of the three-pion exchange, and a repulsive core is introduced 
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by sharply cutting off these parts at the internucleon distance r ~ 0.8 fm. 
At short distances, the interaction is simply given by a constant (but energy- 
dependent) soft core. The Bonn potentials [8] are based on multiple one-boson 
exchange interactions and the two-pion exchange potential calculated in per- 
turbation theory. The two-pion exchange contribution is then approximated 
by an energy-independent u-meson exchange term. Smooth form factors (with 
typical form factor masses ~ 1—2 GeV) cut off the potential at short distances. 
The short range repulsion originates from the u exchange. The second genera- 
tion CD-Bonn potential [9,10] refines the treatment of the one-boson exchange 
amplitudes to take into account the non-local structure arising from the co- 
variant amplitudes. The Nijmegen potentials [11] consist of multiple one-boson 
exchange parts (with separate local and non-local versions), where the inter- 
action parameters of the heavier mesons depend on the partial waves. At very 
short distances, the potentials are regularized by exponential form factors. 
Finally, the local Argonne potential [12] consists of the OPE interaction reg- 
ularized at short distances, and a phenomenological parametrization at short 
and intermediate distances. The core is provided by Woods-Saxon functions 
that are effective at a distance r ~ 0.5 fm. 

Over the past decade, there has also been a great effort in deriving low- 
energy nucleon-nucleon interactions in the framework of effective field the- 
ory (EFT) [13,14,15,16,17,18,19] (for a review see [20]). The EFT approach 
is based on local Lagrangian field theory with low-energy degrees of freedom, 
constrained by the symmetries of QCD. The Lagrangian contains the nucleon 
and pion fields and all possible interactions consistent with chiral symmetry. 
At low energies, the heavy mesons and nucleon resonances can be integrated 
out of the theory. Their effects are encoded in the renormalized pion exchange 
and scale-dependent coupling constants of model-independent contact interac- 
tions. By formulating a power counting scheme, the number of couplings in the 
effective Lagrangian can be truncated and fitted to a set of low-energy data. 
Once the low-energy couplings are determined, the power counting scheme 
enables the EFT to predict other processes with controlled error estimates. 

In addition to the high precision phenomenological interactions, we study the 
Idaho potential [19], which is a chiral model inspired by EFT. The Idaho po- 
tential introduces some model dependence in contrast to the rigorous EFT 
interactions by, e.g., including several terms in the expansion which are sub- 
leading in the power counting with respect to omitted ones. This leads to a 
better description of the experimental phase shifts, comparable in accuracy to 
the realistic models. Since we will argue that our results are dependent on the 
accurate reproduction of the scattering phase shifts, the results for the Idaho 
EFT potential are presented here as well. 

The couplings and parameters of the different potentials described above are 
fitted to the elastic nucleon-nucleon scattering phase shifts over laboratory 
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Fig. 1. The diagonal momentum-space matrix elements of the different high precision 
potentials Vnn versus relative momentum in the 1 So and 3 Si partial wave. 

energies < 350 MeV and the low-energy deuteron properties. As a conse- 
quence, the bare interactions are constrained by the experimental data up to 
a corresponding relative momentum scale of 4 



This is manifest in Fig. 1, where we observe that various high precision poten- 
tial models of Vnn have quite different momentum-space components, despite 
their common treatment of the OPE interaction and the reproduction of the 
same low-energy data. This indicates that the low-energy observables are not 
sensitive to the details or different assumptions of the short-distance physics. 
This insensitivity is a consequence of the separation of long and short- distance 
scales in the nuclear force, and implies that the nucleon-nucleon interaction is 
amenable to an effective theory or renormalization group treatment. 

1.1 The two-body problem in effective theories 

In this work, we propose a renormalization group (RG) approach to the 
nucleon-nucleon interaction which removes the short-distance model depen- 
dencies of the high precision potentials, while preserving their high accuracy 
description of the nucleon-nucleon scattering data and deuteron properties. As 
our approach is a compromise between an EFT treatment and the high preci- 
sion potentials, we first discuss in general terms how effective theory methods 



2.1 fm 



-i 




data 



4 Here, we use conventional scattering units where c = h = H 2 /m = 1, i.e., the 
laboratory energy is related to the relative momentum through E\ a b = 2 k 2 . 



can be applied to the two-nucleon problem. 



The fundamental principle underlying modern effective theory methods is sim- 
ple: the low-energy physics is not sensitive to the details of the high-energy 
dynamics. In other words, one can construct arbitrarily many theories, which 
have the same long- wavelength structure and lead to identical low-energy ob- 
servables, but differ at short distances and higher-energy scales. This ambigu- 
ity can be used constructively. Namely, it is possible to replace the detailed 
short-distance dynamics by simple effective interactions, which preserve the 
low-energy symmetries of the full theory. 

For example, a forerunner to the modern EFT approach is the well-known 
effective range theory of Bethe and Longmire [21]. This allowed a quantitative 
description of low-energy nuclear phenomena without a detailed knowledge 
of the nuclear force. The low-energy scattering amplitude T(k, k; k 2 ) is given 
by two phenomenological low-energy constants independent of an assumed 
nuclear force model [22] 

= --lr k 2 + O(k 4 ), (2) 



T(k,k;k 2 ) a s 2 

where a s denotes the scattering length and r is the effective range parameter. 
The effective range theory encompasses the fundamental principle of modern 
effective theories, as any force models with different short-distance details can 
be tuned to give the same low-energy physics, which is captured in the effective 
range parameters. 

More generally, this insensitivity can be exploited when the short- distance 
dynamics are either poorly understood, or are too complicated for calculation. 
Consider a quantum system described by the Hamiltonian 

H = H + V L + Vk, (3) 

where Vl primarily couples to low-energy states, in our case the long-range 
part of the OPE interaction, and Vu includes the remaining complicated or 
unknown short-distance part of the interaction. 

A first step in building an effective theory is to impose a cutoff A on the inter- 
mediate state energies and momenta. In the modern viewpoint, the cutoff takes 
on a physical meaning. It divides the low-energy states, which are essential 
to the low-energy physics, from the high-energy modes. In this way, one ex- 
plicitly maintains only dynamics that are well-understood. It is, however, not 
possible to simply neglect the remaining Vu, since its effects on the low-energy 
physics need to be included in the form of correction, or so-called countert- 
erms. The correction terms are constrained by demanding that all low-energy 
spectra, low-energy amplitudes, etc., calculated using the bare Hamiltonian H 



5 



are exactly reproduced using an effective Hamiltonian H\ ow k in the truncated 
Hilbert space. 

Schematically, in perturbation theory all amplitudes are of the form 

(Mi) = (f\v hare \i) + ± m^p^ + oiySm)i (4 ) 

n=0 * n 

where the intermediate state summations are over both low and high-energy 
states. We now define an effective low momentum potential V\ ow k by demand- 
ing 

(f\A\i) = (f\V lowk \i) + £ m^M^^ + o( Vl l k ), (5) 

n=0 £/j tl ' n 

where we introduce the correction terms through 

V, owk = V L + 5V ct . (6) 

The physical amplitudes cannot depend on the way in which we split the 
Hilbert space and must therefore be independent of the cutoff. Consequently, 
V\ ow k changes or "runs" with the cutoff in order to cancel the cutoff depen- 
dence arising from the truncated intermediate state summations. This means 
that the effective theory must be RG invariant. Demanding RG invariance 
implies a RG equation for the low momentum interaction 

A</|^)=0 =► ^Mowk = /3([^owk],A), (7) 

where the (5 function of the problem depends on the cutoff-dependent low 
momentum potential, as well as explicit cutoff dependence from the regular- 
ization. Starting from a bare interaction defined in the full Hilbert space, i.e., 
with a large cutoff, the RG equation can be solved to obtain the physically 
equivalent effective potential. This process is referred to as integrating out, 
or decimating the high-energy degrees of freedom. It can be viewed as fil- 
tering out the details of the assumed high-energy dynamics, while implicitly 
incorporating their detail-independent effects on the low-energy physics. 

Since the virtual high-energy states propagate over distances of the order 1/ A 
or shorter, the correction terms are well approximated by contact interactions, 
regardless of the detailed form of V^. The general form of 5V ct should there- 
fore be given by a sum over all local contact operators consistent with the 
symmetries of the problem. In the specific case of the two-body system, we 
have schematically (for an excellent discussion of the details of the approach, 
see Lepage [23].) 



6 



<%(r) = C ± 5 3 (r) + C 2 V 2 5 3 (r) + C 2 ± V • 5 3 (r) V + . . . 

+ ^A^V 2 ^ 3 (r) + ..., (8) 

with dimensionless couplings Cm- To a low-energy probe, the effective theory 
is indistinguishable from the underlying theory, provided the decimation has 
been performed exactly. 

The power counting scheme of the EFT enables one to reliably truncate the 
correction terms, Eq. (8), and fit the couplings directly to a set of low-energy 
data, since a RG decimation of the underlying fundamental theory is at present 
not feasible. In our current approach, we will assume that the nucleon-nucleon 
interaction in the full Hilbert space is given by a high precision potential 
model, and then perform the RG decimation to low momenta exactly. In this 
way, the resulting low momentum potential contains all necessary countert- 
erms to maintain exact RG invariance of the theory, and therefore reproduces 
the experimental phase shift data and the deuteron properties over the same 
kinematic range as the conventional models. 



1.2 A schematic model 



Before considering realistic nucleon-nucleon potential models, we apply the RG 
approach to a simplified schematic model that allows for a straightforward 
solution of the RG equation to construct V\ ow k- In addition, the schematic 
model nicely illustrates the main point of our current approach. Independent of 
the schematic model constructed in the full Hilbert space, the RG decimation 
leads to a model-independent low momentum interaction V\ ow k- 

We consider a separable interaction, which we construct to approximately 
fit the experimental neutron-proton phase shift data in the 1 S channel. For 
this purpose, we found the following separable interaction to be particularly 
useful 5 

V h , IC (k',k)=L(k')R(k) (9) 
= 9 ( 9 ag(k 2 + m 2 ) n \ 

(k 12 + m 2 ) n \ (k 2 + m 2 ) n (k 2 + ri 2 m 2 ) 2n J' 1 ' 

For the diagonal matrix elements, the analogy to the nucleon-nucleon interac- 



5 For simplicity, we have chosen different functions for L(k') and R(k), in order 
to include a repulsive part in the potential without having to introduce several 
separable terms, which would complicate the discussion. 
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I II 

n 1/2 1 

a 1.8199 4.5354 

g [W 1 - 4 ^/ 2 ] 3.1429 3.3542 

m [fm- 1 ] 1.0962 1.3892 

Table 1 

The parameter sets of two separable models of type Eq. (10). 
tion can be seen better 

V h&Ie (k, k) = - (p+m2)2 „ + (A;2 + r?2m 2)2n- ( U ) 

The separable model consists of an attractive part and a short-range, 77 > 1, 
repulsive part. The strength of the repulsion, for given 77, can be constructed 
to fit the 1 S phase shift data reasonably. 

For separable potentials, the Lippmann-Schwinger equation can be resummed 
explicitly and one finds [22] 

L(k')R(k) , x 

2 q° Ljp) Rjp) 
1 + -VJ — 2 — y^-p 2 dp 

7T P — K 

We choose a typical rj — 2, which approximately reflects the ratio of the repul- 
sive to attractive meson contributions m^/2m 7r . The remaining parameters 
of our model, g, m and a, are then fitted to reproduce the effective range 
expansion, 

Tb are (k, k; k 2 ) = — j — , (13) 
^7 ~ 2 r ° ft 

and the change in sign of the experimental phase shifts T bare (/t, n; k 2 ) = at 
1.79 fm" 1 for a given power n. We constructed two exemplary separable 
models with the parameters given in Table 1. The separable models are fitted 
to a neutron-proton 1 S scattering length, a s = — 23.73 fm, and effective range 
r = 2.70 fm. 

In Fig. 2, we show the phase shifts calculated from the models of Vb are . We 
observe that the separable models I and II can achieve a reasonably realistic 
description of the nucleon-nucleon phase shifts for laboratory energies below 
-Slab < 350 MeV, i.e., for relative momenta k < 2.1 far 1 . 

As discussed in the previous Section, the low momentum amplitudes can be 
preserved with an effective interaction Vj ow k acting in the Hilbert space of 
low momentum modes, k < A, exclusively. The low momentum potential is 
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Fig. 2. The left figure shows the phase shifts calculated from the separable mod- 
els I and II in comparison with the effective range expansion and the Nijmegen 
multi-energy phase shift analysis (PSA93) [24,25]. In the right figure, the phase 
shifts calculated from the corresponding low momentum interactions are given. The 
low momentum interactions are obtained by solving the RG equation, Eq. (14), to 
a scale A = 2.1 fm -1 . 

renormalized by the high momentum modes according to an RG equation, 
of general form given by Eq. (7). In the current approach, we demand that 
the low momentum half-on-shell T matrix (and consequently the phase shifts) 
are RG invariant, see also [3,4]. The RG equation which follows from this 
requirement is derived in [26]. For separable interactions, the RG equation for 
the diagonal components of the low momentum interaction simplifies and is 
given in closed form by 6 



d 2/tt V lawk (k,k;A)V ]owk (A,A;A) 

— V low k {k, k; A) - - _ (fc/A)2 



1 + ^A1^«A) 

7T P — A A 



(14) 



The initial condition of the RG equation is given by the bare interaction, 



Vlow k(fc, fcl A ) = Vbare^, k), 



(15) 



at an initial scale A where the high momentum modes are negligible, 



VLare(Ao,A ) ~ 0. 



(16) 



6 The RG equation is regularized by imposing a small gap between the low momen- 
tum Hilbert space, k ^ A — e, and the cutoff A in the integration over intermediate 
momenta, p\ 00 p ^ A. In this way the T matrix can be calculated to chosen accuracy. 
In other words, the effect of this regularization on the scattering amplitude is O(e). 
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Fig. 3. The diagonal momentum-space matrix elements of Vbare and Vi ow k derived 
from the separable models I and II. As a comparison, we show the results for the 
diagonal matrix elements of the V[ow k from Section 3 derived from the high precision 
nucleon-nucleon interactions. 

The RG equation is then used to evolve the bare potentials to low momenta, 
thus removing the model-dependent realization of the high momentum parts. 
We solve the separable RG equation to a cutoff A = 2.1 fm -1 , which corre- 
sponds to a laboratory energy scale up to which the realistic nucleon-nucleon 
interactions are constrained by experiment, ~ 350 MeV. First, we com- 
pare the phase shifts calculated in the low momentum space 



tan S = —k T low k (k, k; k ) 
2 , V\ owk (k,k;A) 



2 ^ .^A) 
7r o V ~ k 



(17) 
(18) 



and verify in Fig. 2 that the phase shifts are indeed reproduced by Vi ow k- In 
the RG approach, we have by construction 



^low k(k , k\ k ) — Tbare(^ ; k] k 



(19) 



In Fig. 3, we show the resulting diagonal matrix elements of V\ ow k in compari- 
son to the V\ ow k derived from the high precision nucleon-nucleon interactions. 
Although the bare separable models I and II differ substantially, we find that 
the low momentum interactions in the truncated Hilbert space are identical. 
Furthermore, we observe that, to a surprisingly good accuracy, these simple 
models can in fact reproduce the diagonal matrix elements of the low mo- 
mentum interaction obtained from the realistic potential models. In the RG 
equation for separable models, the matrix element corresponding to zero phase 
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shift remains invariant, 

Mow k(«, K| A) = Vbarc(/«, «) = T(k, K\ K 2 ) = 0. (20) 

Interestingly, we find that the diagonal matrix element of the realistic V\ ow k 
also changes sign approximately at relative momentum k. 

1.3 Organization 

This work is organized as follows. In the next Section we discuss the calcu- 
lational methods used for constructing the effective low momentum interac- 
tion Vjowk- in particular, we show how Mow k is related to the solution of the 
Bloch-Horowitz equation. Since the equivalence of the model space effective 
interaction methods and the RG approach is discussed in detail in [26], we will 
focus on the discussion of the effective interaction methods used in this work. 
In Section 3, we present the results for the low momentum interactions derived 
from the realistic potential models. We give results for all partial waves with 
J ^ 4, which provides ample evidence for the claimed model independence of 
Mow k- By studying specific partial waves, we will show how the model inde- 
pence is reached as the cutoff is varied, and also demonstrate its relation to 
the phase shift equivalence and the common long-range pion physics of the dif- 
ferent potential models. Finally, we verify that the low momentum interaction 
indeed reproduces the scattering data below _E lab < 350 MeV and the deuteron 
binding energy with similar accuracy as the high precision potentials. After 
summarizing our results, we elaborate on the advantages of using Mow k as the 
microscopic input interaction for nuclear many-body calculations. 



2 Calculational methods 

The schematic model provides a nice example illustrating the use of nonpertur- 
bative RG techniques to construct the low momentum effective theory. One 
could therefore apply the same direct method as for the schematic model. 
However, in the general non-separable case, the direct solution of the RG 
equation would be more involved. Therefore, we have adapted conventional 
effective interaction techniques, also referred to as model space methods, for 
the calculation of the low momentum interaction in the general case. These 
methods have been successfully used to derive nuclear shell model interactions 
in a truncated Hilbert space, see e.g., the review [27]. Both the RG and model 
space techniques are concerned with integrating out the high-energy states, 
such that the low-energy observables remain invariant in the simpler effective 
theory. In spite of their similarities, the formal equivalence of conventional 
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effective interaction theory and the RG was only recently demonstrated for 
the two-body problem [26]. There, we have shown that the low momentum 
Hamiltonian obtained from the solution of the RG equation is equivalent to the 
effective theory derived using Bloch-Horowitz or Lee-Suzuki projection meth- 
ods. Since our current work makes use of this equivalence for the calculation 
of Viowk, we briefly discuss the relevant techniques below. 



2.1 Integrating out the high momentum modes by similarity transformations 



Consider a physical system, for which we are only interested in its low-energy 
properties. The first step in effective interaction theory is to define Feshbach 
projection operators onto the physically important low-energy model space, 
the so-called P-space, and the high-energy complement, the Q-space, 



^ = EI0(<I= ^d Q= £ m\=\ |. (21) 

i=\ \ / i=d+l 

The projection operators satisfy P + Q — 1, PQ = QP = 0, P 2 = P and 
Q 2 = Q, and we have written the projection operators in matrix form as well. 
In most applications, the eigenstates of the unperturbed Hamiltonian H are 
taken as basis states 

The Schrodinger equation can then be written in block form as 





PHP 



phq\ 



P\V)\ P\V) 

' E\ I . (22) 



[qhp qhq) \Q\*)) \Q\V) 



By using the second block-row of Eq. (22), the Q-space projection Q can 
be eliminated and one obtains the fundamental equations of Bloch-Horowitz 
effective interaction theory [28,29] 

PHZ k (E)P\V) = EP\V), (23) 

where the effective Hamiltonian TC™ k (E) is obtained by the solution of the 
Bloch-Horowitz equation 

H ™^ = H + H E^QHQ H - < 24 > 

The Bloch-Horowitz equation generates an effective Hamiltonian, whose P- 
space projection H^ k (E) = PW™ k(-^0 P ls operative only within the low- 
energy model space and, e.g., enables an exact diagonalization in shell model 
applications. The Q-space states have been decoupled from the problem in 
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a way that preserves the low-energy spectrum exactly, i.e., the d lowest-lying 
eigenvalues. The eigenstates of H^ k (E) are simply given by the P-space pro- 
jections of the exact eigenstates P |\P n ). The effective Hamiltonian H^ k (E) 
depends on the exact energy eigenvalue one is solving for, and thus necessi- 
tates a self-consistent treatment. Although the self-consistency problem can 
be solved nicely using the Lanczos algorithm [30], the use of energy-dependent 
two-body interactions in many-body calculations can pose computational and 
conceptual difficulties. 

A refinement of the Bloch-Horowitz theory that avoids these difficulties is 
the Lee-Suzuki similarity transformation [31,32]. The Lee-Suzuki method con- 
structs a similarity transformation that brings the Hamiltonian to the follow- 
ing block structure 



& 1 HQ = H{f w k = 



PHP PHQ" 
QHQj 



(25) 



Since the determinant of a block-triangular matrix factorizes into the deter- 
minants of the independent P-space and Q-space blocks, the low-energy spec- 
trum can be calculated exactly from the smaller effective model space Hamilto- 
nian k = PH\^ k P. As in the Bloch-Horowitz approach, the eigenstates 
of -ffiowk are given by P|\I/ n ). Making an ansatz for the Lee-Suzuki similar- 
ity transformation, one parametrizes O as a non-orthogonal transformation in 
terms of the so-called wave operator uj 



Q = 1 + uj 



and 1 = 1 — uj 



' 1 0^ 



(26) 



where uj = Q uj P connects the Q and P-space. Inserting the ansatz into 
Eq. (25) leads to a non-linear constraint on uj, the decoupling equation, 

QW^ k P=0 =>• QHP + QHQuj-ujPHP-ujPHQuj = 0. (27) 

Once a solution to this non-linear operator equation is found, the energy- 
independent effective Hamiltonian is given by 

< k = PHP + PHQujP = PHP + PVQujP, (28) 

where we have assumed that P and Q commute with the unperturbed Hamil- 
tonian. 

The relation to the Bloch-Horowitz theory can be clarified by inserting the 
Lee-Suzuki solution, Eq. (28), into the decoupling equation, Eq. (27), and 
solving for the wave operator. The resulting formal expression for the wave 
operator allows one to express the energy-independent Lee-Suzuki k as a 
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summation of the lowest d, self-consistent Bloch-Horowitz solutions H™ k (E), 
where each term is weighted by the projection of the corresponding eigenstate 

of the effective theory P \V n ) ((* n | p) 7 

*S > = t («"> + ™ ^qWq Hp) P |*„) (MP) 

= •£ H™ „(£„) P | *„>((*„| P). (29) 

n=l 

In this work, we apply the Lee-Suzuki formalism to the free-space nucleon- 
nucleon problem to derive the energy-independent low momentum interaction 

View k = #iow k - PH P = PVP + PVQ uo P. (30) 

This approach has been proven to be equivalent to the solution of the RG 
equation that results from imposing a momentum-space cutoff and demanding 
RG invariant half-on-shell (HOS) T matrices [26]. For a given partial wave, 
the P and Q operators are defined in a continuous plane wave basis as 

2 pA 2 r°° 

P = - p 2 dp\p)(p\ and Q = - q 2 dq\q)(q\. (31) 

IT JO IT J A 



There are several methods to solve the non-linear decoupling equation for the 
wave operator. For the two-body problem, e.g., the exact solutions for the 
eigenstates can be used, and Vj wk can be calculated directly using Eq. (29). 
For more complex applications, one solves the decoupling equation using iter- 
ative techniques. Since the iterative methods are numerically very robust, we 
have utilized the iterative algorithm of Andreozzi [33] to calculate V\ ow k- The 
defining equations of this iteration scheme are given by 

Xo = -(QHQ)~ 1 QHP (32) 

Xn = (QK-i)) _1 # n _i k (w n _i), for n > 1, (33) 

where 

u n = X + X 1 + . . . + X n (34) 

Q{uj n ) = QHQ - Q u n PHQ (35) 

and 

#iow kK) = PHP + PHQ u n P. (36) 

7 The tilde denotes the bi-orthogonal complement, since projections of orthogonal 
\^ n ) are not necessarily orthogonal. 
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The iteration converges when H^ k (uj n ) xs H^ k (u; n -i). 

A convenient feature of the algorithm is that each iteration only requires one 
matrix inversion in the Q-space. The Q-space matrix inversion is manageable, 
typically requiring a set of ~ 60 Gauss-Legendre mesh points extending to 
25 fm" 1 to include the high momentum behavior of the input interaction. 
Finally, we mention that variations of the algorithm exist in which all matrix 
inversions occur in P-space, which can offer significant computational advan- 
tages for large-scale problems [33] . 

2.2 A refinement: hermitian low momentum interactions 

By construction, the Lee-Suzuki k is non-hermitian because the eigen- 
states of the effective theory are projections of the full eigenstates onto the low 
momentum subspace. The non-hermiticity is a result of eliminating the energy- 
dependence of the Bloch-Horowitz theory by means of the non-orthogonal 
similarity transformation. 8 However, one can perform an additional similar- 
ity transformation on k to obtain an energy-independent and hermitian 
Hamiltonian [34,35], which we denote by H low k and correspondingly, the her- 
mitian low momentum interaction by V\ ow k . Both similarity transformations 
combine to give a unitary transformation, which is equivalent to the well- 
known Okubo transformation [36]. 

Here, we follow the method of Andreozzi [33] in constructing the second sim- 
ilarity transformation. Using the definition of the Lee-Suzuki transformation, 
Eq. (25), one finds 

e^^^K^e 1 ©, (37) 

where we have used H = W. Projecting onto the low momentum space, the 
P-space block can be written as 

(P + J H^ k = (H^ k )\p + J u). (38) 

The operator P + or a; is hermitian and positive definite. Therefore, it can 
be written by means of a Cholesky decomposition as a product of a lower- 
triangular matrix L and its adjoint LL Writing 

(P + Juj) = LL ] (39) 

in Eq. (38) and multiplication by L^ 1 from the left and from the right, 

8 Although the Bloch-Horowitz Hamiltonian H^ k (E n ) is hermitian for a given 
energy E n , one encounters the same non-orthogonality problem between eigenstates, 
since the Hamiltonian changes self-consistently depending on the energy E n . 
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one finds 



L' 1 [n\l^L = Lt^ k (Lt)- = (L- (h\1^l)\ (40) 

Therefore, we can write the hermitian h\^ w k directly in terms of the Lee- 
Suzuki transformation and the Cholesky decomposition as 

^lot^^^lowk^)- 1 , (41) 

and the combined transformation U = {L^)^ 1 is unitary. 

We mention that the differences between the non-hermitian V\ ow k and the 
hermitian V\ ow k can be understood in the RG approach. The non-hermitian 
effective theory is obtained by introducing a cutoff in momentum space and 
imposing RG invariance for the HOS T matrix with an energy-independent 
low momentum potential. Since the initial conditions of the RG equation are 
given by the bare interaction, this implies that the effective theory preserves 
the low momentum components of the full wave functions, in addition to 
the phase shifts and the bound state poles. The choice of HOS T matrix 
invariance is motivated by the physical observation that the low momentum 
components of the wave functions mainly probe the long-range OPE part of 
the interaction. Conversely, the hermitian effective theory results from RG 
invariant on-shell T matrices. Thus, the resulting hermitian V\ ow k no longer 
preserves the low momentum components of the full wave functions, although 
the observable phase shifts and bound state poles are preserved. In practice, 
we will find that the numerical differences between the hermitian and non- 
hermitian interactions are quite small. 

Finally, we note that the Lee-Suzuki similarity transformation method is an 
algebraic realization of the more general folded diagram method of Kuo, Lee 
and Ratcliff [27,37,38]. In this connection, one can recast the results of this Sec- 
tion in a diagrammatic formalism, which can be useful in many-body systems. 
Moreover, the folded diagram formulation has been applied to the two-body 
problem to prove diagrammatically that the low momentum effective theory 
preserves the low-energy scattering observables of the input interaction [3] . 



3 Results 

We now discuss the main results for the low momentum interaction Vi owk , 
which we have obtained by integrating out the high momentum components 
of different realistic potential models. Most of the results shown here are for 
a momentum cutoff of A = 2.1 fm -1 , which corresponds to a distance scale 
1/A ~ 0.5 fm. With cutoffs chosen in this range, the model-dependent pa- 
rameterization of the repulsive core at r < 0.5 fm in the bare interactions is 
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not resolved in the effective theory. Consequently, one is able to reproduce 
the deuteron pole and scattering observables below E\ a ^ = 350 MeV without 
introducing a specific short-distance model. Moreover, the RG invariance of 
the decimation guarantees that our conclusions are independent of the precise 
value of the cutoff. 

Before we present our results in detail, we summarize our main findings and 
refer the reader to the key figures. The most important properties of V\ ow k 
are: 

1- Mow k renormalizes to a nearly universal interaction as the cutoff is lowered 
to the scale of A ~ 2.1 fm _1 . We will argue that the model-independence 
of the diagonal matrix elements of V\ ow k, as shown in Fig. 4 and Fig. 5, 
is driven by the phase shift equivalence of the input models. The cor- 
relation between the phase shift equivalence and the collapse of the low 
momentum interactions is clearly demonstrated in Fig. 8. We also find 
practically identical results for the hermitian Ki ow k, see Fig. A.l and 
Fig. A.2. 

2. V\ ow k reproduces the phase shifts and the deuteron binding energy with 
similar accuracy as the high precision models. The phase shifts calculated 
from V\ ow k are shown in Fig. 16 and Fig. 17. 

3. In the effective theory one does not have to make model-dependent as- 
sumptions on the short-distance interaction. The RG evolution filters 
out the short- distance details of the input potentials, while preserving 
the model-independent effects of the high momentum modes on the low 
momentum observables, and leaving the long-range part (OPE) of the 
interaction intact. This can be seen most clearly in Fig. 11, Fig. 12 and 
Fig. 13. 

We also point out that the numerical differences between V\ ow k and the hermi- 
tian V\ ow k are quite small, see Fig. 15. For brevity, we mostly present results 
for V\ ow k in this work, since the V\ ow k is easily derived from V\ ow k using the 
transformation described in Section 2.2. 



3. 1 Model-independence of the low momentum interaction 

For cutoffs around A ~ 2.1fm _1 , the effective theory preserves the nucleon- 
nucleon phase shifts up to laboratory energies of E^b ~ 350 MeV. This corre- 
sponds to the scale of the scattering data that all the high precision potential 
models are fitted to. We have argued that, based on the separation of mass 
scales, the model-dependence of the realistic interactions should be reduced as 
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Fig. 4. Diagonal momentum-space matrix elements of the Vj ow k obtained from the 
different potential models for a cutoff A = 2.1 fm -1 . Results are shown for the 
partial waves J ^ 4. 
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Fig. 5. The diagonal matrix elements of V\ ow k are continued from Fig. 4. 



we integerate out the ambiguous high momentum modes. In Fig. 4 and Fig. 5, 
we present the results of the RG decimation. We observe that the diagonal 
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Fig. 6. The off-diagonal momentum space matrix elements of the Vi ow k obtained 
for a cutoff A = 2.1fm _1 (left) and A = l.Ofm - (right) in the 1 So partial wave. 

matrix elements of the low momentum interactions derived from the different 
potential models are nearly identical, and thus insensitive to the input poten- 
tial model for cutoffs A < 2.1 fm -1 . The largest relative differences between 
the Viowk are found in the 3 F 2 , 3 D 3 , ^3 and 3 F 3 partial waves. These differ- 
ences are correlated with the deviations in the phase shifts, i.e., the different 
accuracies, of the potential models. We will discuss this point in detail below. 



We note that if we perform the additional similarity transformation to obtain 
the resulting hermitian V^iow k, we find the same model-independence and prac- 
tically identical diagonal matrix elements, as shown in Fig. A.l and Fig. A. 2 in 
the Appendix. Since the hermitian V]^ k is obtained by symmetrizing V[ ow k, it 
is expected that the effects on the diagonal matrix elements are not significant. 

One obtains further insight into the collapse of V\ ow k by studying the off- 
diagonal matrix elements. These are shown in the 1 So partial wave in Fig. 6. 
We find a similar nearly universal behavior for the off-diagonal matrix ele- 
ments of V[ow k, although those calculated from the Idaho potential differ at 
k' ~ 1.2 fm _1 , which is approximately the mass scale given by the 2ir exchange. 
The conventional models include 2tt contributions in which non-nucleonic in- 
termediate states are explicitly taken into account. In contrast, the Idaho po- 
tential starts from a local, effective Lagrangian, in which nucleons and pions 
are the only explicit degrees of freedom. For example, contributions from A 
intermediate states between two successive pion exchanges are contracted to 
a renormalized C^-^ vertex. Moreover, the conventional potential models all 
contain a strong, short-range repulsion associated with the exchange of heavy 
mesons, mostly originating from the u exchange. In the Idaho potential, the 
repulsive core is not explicitly present. However, the effects of the short-range 
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Fig. 7. Comparison of the HOS T matrices calculated from the different potential 
models in the 1 So partial wave. 

repulsion on low-energy observables are contained in contact terms that are 
fitted to the low-energy data [19]. Therefore, in momentum space, the Idaho 
potential changes over to contact terms at momentum transfers above the 
range given by 2m n . 



While the Idaho potential is phase shift equivalent to the conventional mod- 
els, the HOS T matrices, i.e., the low momentum components of the wave 
functions, differ significantly from the conventional models at energies and 
momenta larger than the above scale of 2m n . In Fig. 7, we find that the con- 
ventional models give remarkably similar HOS T matrices over their range of 
phase shift equivalence. This approximate HOS T matrix equivalence can be 
understood by expanding the HOS T matrix around the fully on-shell part. 
The difference can be written as a wound integral that contains the off-shell 
behavior, 

k 2 — k! 2 f°° / \ 

T(k', k; k 2 ) = T(k, k; k 2 ) + - J dr sm(k'r) (u k (r) - v k (r)) , (42) 

where Wfc(r) = r\l/ fc (r) denotes the exact and v k (r) is the asymptotic S-wave 
wave function 

lim Mfe(r) — > vdr) ~ sin(A;r + 5). (43) 

r— >oo 

We can use Eq. (42) to express the difference of the HOS T matrices for two 
different, but phase shift equivalent bare interactions. This leads to 

k 2 - h' 2 r R 

ATM, k; k 2 ) = / dr sin(JfeV) Att fc (r), (44) 

k' Jo 

where R is the distance over which the two models differ substantially at 
short distances, which we can take for simplicity as the size of the repulsive 
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Fig. 8. Comparison of the phase shifts given by the Vnn models and the diagonal 
matrix elements of V\ ow k in the 3 D3 and 3 F3 partial waves. Results for Vj ow k are 
shown with a cutoff A = 2.1 fm _1 . 

core. Regardless of the details of the differences, at scattering energies well 
below the strength of the short-range repulsion, the differences between the 
wave functions in the core are suppressed by strong two-body correlations. 
Therefore, it is physically reasonable that the conventional realistic potentials 
give similar HOS T matrices at low energies. 

The deviations in the off-diagonal matrix elements of V\ ow k derived from the 
Idaho potential are obviously removed, if we integrate out further below 2m 7r . 
As shown in Fig. 6, we find that the off-diagonal matrix elements in the 1 So 
partial wave collapse onto one curve as well for A < 1.2 fm~ . 

The preceding observations strongly suggest that the collapse of the diagonal 
matrix elements of V\ ow k is driven by the phase shift equivalence of the input 
models, while the collapse of the off-diagonal matrix elements is controlled by 
approximate HOS T matrix equivalence at low energy and momentum. The 
latter can in turn be mainly attributed to a common off-shell behaviour in the 
OPE interaction, with some differences arising from the two-pion exchange and 
the crossover to contact interactions. The correlation between the diagonal 
matrix elements and the phase shifts is nicely shown in the 3 D 3 and 3 F 3 
partial waves, where we observe the largest relative deviations from model- 
independence. The left panels of Fig. 8 clearly show that the older generation 
Paris and Bonn potential models give significantly different phase shifts in 
these partial waves. This leads to the observed differences in the Vj ow k derived 
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from these bare interactions, as can be seen in the right panels of Fig. 8. 
For comparison, we have also shown the results of the Nijmegen multi-energy 
phase shift analysis (PSA93) [24,25]. 



A further example of this correlation can be seen in Fig. 9, where the Idaho 
potential leads to considerably different phase shifts in the 3 F 2 partial wave 
above laboratory energies > 100 MeV. This corresponds to a relative mo- 
mentum of k > l.Ofm - . Accordingly in Fig. 10, we find that the diagonal 
matrix elements of Vj ow k derived from the Iadaho potential deviate from the 
others for momenta larger than this scale. As for the off-diagonal matrix ele- 
ments in the 1 So partial wave, these differences vanish as the cutoff is lowered 
to A < l.Ofm -1 . 
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Fig. 11. The collapse of the diagonal momentum-space matrix elements of Vj ow k as 
the cutoff is lowered to A = 2.1 fm -1 in the 1 So partial wave. 

In the N 2 LO version of the Idaho potential used here, the contact interactions 
do not contribute to the L ^ 3 partial waves. Thus, in the 3 F2 partial wave the 
Idaho potential is completely given by pion exchange. Our results demonstrate 
that the effects of the short- distance physics cannot simply be ignored. Rather, 
their renormalization effects on the low momentum physics has to be taken 
into account. If we lower the cutoff to A < 1.0 fm -1 , the effects of the high 
momentum modes is sufficiently small, and one finds that indeed the physics 
is given adequately by pion exchange. This can be understood in the Born 
approximation, which is a reasonable approximation for the high partial waves, 
especially at low energies. In the Born approximation, the integral term of 
the Lippmann-Schinger equation is neglected, and thus the low momentum 
interaction in the cutoff Hilbert space is given by simply cutting off the bare 
interaction. 



The RG approach has the advantage that the nature of the low momentum 
interaction can be revealed by studying the change of V[ ow k with the cutoff 
scale. In this way, the RG evolution clearly provides evidence for the separation 
of scales in the two-nucleon problem, see also [3]. Furthermore, we will show 
that it also provides physical insight into the low momentum interaction. 

In the evolution of the 1 So matrix elements shown in Fig. 11, we find that the 
model independence of V\ ow k is reached at somewhat larger values of the cutoff 
around A ~ 3.5 fm -1 . This scale lies below the u and p meson masses, and 
consequently the details of the repulsive core are not resolved in the effective 
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Fig. 12. The collapse of the diagonal momentum-space matrix elements of V\ ow k as 
the cutoff is lowered to A = 2.1 fm -1 in the 3 Si partial wave. 
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Fig. 13. The collapse of the diagonal momentum-space matrix elements of Vj ow k as 
the cutoff is lowered to A = 2.1 fm -1 in the coupled 3 Si- 3 Di block. 



theory. Also note that for A > 3.0 fm~ , the Idaho potential is not renormal- 
ized, since it is an EFT potential without high momentum components by 
construction. 
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E D [MeV] Paris Bonn A Nijm. I Nijm. II Arg. u 18 CD Bonn Idaho A 



Vnn -2.2218 -2.2242 -2.2246 -2.2242 -2.2247 -2.2238 -2.2242 
quoted -2.2249 -2.22452 -2.224575 -2.224575 -2.224575 -2.224575 -2.2242 
Vi ow k -2.2218 -2.2242 -2.2246 -2.2242 -2.2247 -2.2238 -2.2242 

Table 2 

Comparison of the deuteron binding energy Ed for V[ ow k and the bare interactions. 
We note that all V\ ow k are derived with the same momentum mesh. For the values 
reported here, we have chosen not to optimize the mesh for each bare interaction 
separately. 
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Fig. 14. Comparison of the momentum-space S and D-state deuteron wave functions 
calculated from Vi ow k> the hermitian Vi ow k and the bare Vnn- Results are shown 
for a cutoff of A = 2.1 fm -1 and the Argonne v%s potential. 

Conversely, we find in Fig. 12 that the collapse in the 3 Si partial wave occurs 
at a comparatively lower value for the cutoff of A ~ 2.1 — 2.6 fm _ . This arises 
from the fact that the tensor force is operative in the triplet S-wave. Since 
Mow k acquires a large second-order renormalization from the tensor force, one 
indeed expects that V\ ow k becomes model-independent at a lower cutoff in the 
3 Si partial wave. 



The nature of the tensor interaction of V[ ow k can be studied by considering 
the couping between the 3 Si and 3 D! partial wave. In this channel, only the 
tensor force enters. The results for the RG evolution in the 3 Si- 3 D! block 
are presented in Fig. 13. We observe that the realistic models are practically 
identical for momenta below k < 0.7 fm -1 , which corresponds to the pion 
mass. This clearly demonstrates that the tensor part of the input models differs 
significantly for momenta above m n . However, the effects of the different high 
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Fig. 15. Comparison of the diagonal (left) and off-diagonal (right) matrix elements 
of V\ ow k and the hermitian V\ ow k in the 1 So partial wave. Results are shown for a 
cutoff of A = 2.1 fm -1 and the Argonne vis potential. 

momentum tensor components on the low momentum physics can be taken 
into account by employing a renormalized interaction, which in turn becomes 
model-independent at a cutoff of A ~ 2.1 fm -1 . 

The RG decimation to V\ ow k preserves the deuteron binding energies of the 
bare interactions. In Table 2, we give the binding energies obtained from the 
V[ow k, and verify that the deuteron poles are in fact reproduced in the effective 
theory. We note that we have not optimized the momentum mesh for each 
input Vnn, which would further improve the agreement with the quoted values. 
Referring to Fig. 14, the low momentum components of the bare deuteron S 
and D-state wave functions are RG invariant as well. Moreover, although the 
hermitian V\ ow k does not exactly preserve the low momentum projections 
of the wave functions, we find in Fig. 14 that the deuteron wave functions 
calculated from V\ ow k and the hermitian V\ ow k are practically identical. 

Finally, we note that the same general results hold for both V\ ow k and V\ ow k, 
and that the numerical differences are quite small. This is demonstrated as 
an example for the l So partial wave in Fig. 15. It can also be seen that the 
differences are as expected most noticeable for the far off-shell parts of the 
low momentum interaction. 

3.2 Phase shift equivalence of the low momentum interaction 

The mo del- independent low momentum interaction reproduces the experimen- 
tal phase shift data with similar accuracy as the high precision nucleon-nucleon 



26 



60 
40 
20 


-20 
4 

r oB 3 
U 

£ 2 

2 1 

3 

w 

1U 



10 



V 150 

\^ 100 




\ " 15 


: 


100 200 300 ^ a 

f*^^ -20 
F E i -30 


100 200 300 1Q 

«k 
^^^e -10 
'P ** * 

1 ^2TS 


« 100 200 300 


100 200 300 

io 

5 


100 200 300 ^ 

...» — * — ^5 
* 3 P 2 0.0 


100 200 -O00 


t 100 200 300 
\ 10 

V S 


100 200 300 


100 200 300 

PSA93 

* CD Bonn 
o hermitian V 



100 

E, 



200 300 

[MeV] 



Fig. 16. The neutron-proton phase shifts and mixing parameters calculated from 
the hermitian V\ ow k in the partial waves J ^ 4 are shown. The V\ ow k is derived 
from the CD Bonn potential with a cutoff of A = 2.1 fm -1 . In comparison, the phase 
shifts and mixing parameters calculated from the bare CD Bonn potential as well 
as the results of the Nijmegen multi-energy phase shift analysis [24,25] are included. 
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Fig. 17. The neutron-proton phase shifts and mixing parameters are continued from 
Fig. 16, for details we refer to the text of Fig. 16. 



potentials. Therefore, V\ ow k fits the measured phase shifts for laboratory en- 
ergies below -Ei a b = 350 MeV with a similar x 2 /datum « 1 — 2, depending on 
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which conventional potential model Vnn is used. Since the hermitian V\ ow k 
preserves this content as well as the hermiticity of the Hamiltonian, we demon- 
strate the high accuracy of the hermitian V\ ow k here. In Fig. 16 and Fig. 17, 
the neutron-proton phase shifts and the mixing parameters in coupled chan- 
nels for all partial waves J < 4 are given. For clarity of the Figures, we show 
these for the low momentum V\ ow k derived from the CD Bonn potential only. 
Since both the hermitian and non-hermitian low momentum interactions pre- 
serve the phase shifts of the input models, results of the same accuracy are 
derived from the various other models. For comparison, the phase shifts and 
mixing parameters calculated from the bare CD Bonn potential and the re- 
sults of the Nijmegen multi-energy phase shift analysis [24,25] are included in 
Fig. 16 and Fig. 17. 

We conclude that the low-energy scattering data can be equally well repro- 
duced by an effective low momentum interaction Vy ow k restricted to the Hilbert 
space of low momentum modes, k < A^ata- This is achieved with a model- 
independent interaction, without introducing further assumptions on the de- 
tails of the short-distance dynamics. 



4 Summary and advantages 

We have shown that the RG decimation of different high precision nucleon- 
nucleon interactions to low momenta leads to a unique interaction, called 
Viow k- This result is obtained by truncating the full Hilbert space to an ef- 
fective space, which consists of momentum components up to the scale of 
the constraining low-energy data. The differences in the realistic interactions 
arise from the assumed high momentum dynamics. By restricting the inter- 
action to the low momentum space, the details of the short- distance physics 
are not resolved, and the detail-independent effects of the high momentum 
modes on the low momentum observables are included in the renormalized 
effective interaction. The existence of a nearly universal V\ ow k demonstrates 
that it is possible to separate the physics of low momentum nucleons from 
the ambiguous short- distance parametrizations used in the realistic potential 
models. 

The RG invariance of the effective theory guarantees that the HOS T matrix 
and the deuteron properties of the full theory are exactly preserved in the 
model space. Consequently, V\ ow k is phase shift equivalent to the high preci- 
sion potential models and reproduces the experimental elastic scattering data 
with similar accuracy. A further refinement transforms V\ ow k into a hermi- 
tian low momentum interaction V\ ow k, which remains model-independent and 
preserves the low-energy on-shell amplitudes of the bare interactions. In this 
manner, the resulting RG decimation incorporates the model-independent ef- 
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fects of the short- distance physics in the low momentum interaction, while the 
model-dependent effects are filtered out. 



We believe that the low momentum interaction Vj ow k satisfies the first of the 
two requirements of a microscopic many-body theory given in the Introduc- 
tion. Namely, that the input interaction is, as best as possible, based on the 
available low-energy data. A necessary condition for the collapse of the low 
momentum interactions is the separation of long from short-distance scales 
in the nuclear problem. In addition, we have argued that the accuracy of the 
reproduction of the phase shifts drives the collapse of the on-shell components 
of the low momentum interactions. Furthermore, the evolution of Vi ow k with 
the cutoff provides physical insight into the separation of scales in the nuclear 
problem, the effects of the short-range repulsion in the renormalization, and 
the tensor content of the low momentum interaction. 



^iow k sums high momentum ladders in free space. Therefore, the low momen- 
tum interaction is considerably softer than the bare interactions and does not 
have a repulsive core. As a consequence, when used as the microscopic input in 
the many-body problem, V\ ow k does not lead to strong high momentum scat- 
tering effects in the particle-particle channel, which necessitates a Brueckner 
resummation or short-range correlation methods for the conventional models. 



The low momentum interaction should be regarded as a new input potential 
for many-body calculations. However, we emphasize that in contrast to the 
Brueckner G matrix, V\ ow k does not have high momentum components above 
the cutoff and therefore one does not have to be concerned about double count- 
ing the high momentum contributions. In EFT, this is understood easily, since 
the ultraviolet divergences can be renormalized in free space and the many- 
body dynamics does not lead to new ultraviolet divergences. More precisely, 
in any many-body diagram, the excitations to high momentum states that 
are absorbed in Vj ow k vanish due to the momentum space cutoff by construc- 
tion. Using V\ ow k is simply analogous to using the relatively soft CD Bonn 
interaction, instead of the stronger Argonne t>i 8 potential. 



Since V\ ow k accounts for a large part of the phase space in the many-body sys- 
tem, it will be interesting to explore whether calculations starting with V\ ow k 
could be organized in a more perturbative fashion. Encouraging results for a 
perturbative expansion of the shell model effective Hamiltonian in few-body 
systems has been demonstrated for the deuteron in [39]. For larger systems, 
this idea more relies on a geometrical phase space argument. 
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A Diagonal matrix elements of the hermitian V\ OVf k 
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Fig. A.l. Diagonal momentum-space matrix elements of the hermitian V\ ow k ob- 
tained from the different potential models for a cutoff A = 2.1 fm -1 . Results are 
shown for the partial waves J ^ 4. 
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Fig. A. 2. The diagonal elements of the hermitian V\ ow ^ are continued from Fig. A.l. 
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